UNCLASSIFIED A K 

security  classification  or  thi»  pace  (m >«»  Dot,  hum  tltS 

REPORT  DOCUMENTATION  PAGE 

/ ...  ~ IToovt  accession  wo. 

' ^RCSL-TR-78^21 

4.  TITLE  (mtd  Subtltlo)  ™” 

I COMPUTER  PROGRAM  FOR  ANALYZING  PROJECTII 


O £ 


I -*  A LUMrUTbR  PROGRAM  F 
/ , , FLIGHT  TRAJECTORY  » 


7. 

I John  MiFerriter  l 
1 JohnJ.lReba  1 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

*■  RECIPIENT'S  CATALOG  NUMBER 


i-f  lYPEOFREPOR'l  W“P*RK>Q~COUABEJ 

1 1 Te(*j"ig^y 

t.  PERFORMING  0*0.  REPORT  NUMBER 
• ■ CONTRACT  OR  ORAHT  NUMBERf*) 


1-  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

/ Director,  Chemical  Systems  Laboratory 
Attn:  DRDAR-CLN-D 
Aberdeen  Proving  Ground,  Maryland  21010 

11.  CONTROLLING  OFFICE  NAME  AND  ADORES* 

Director,  Chemical  Systems  Laboratory 

Attn:  DRDAR-CLJ-R 

Aberdeen  Proving  Ground,  Maryland  21010 

14.  MONITORING  AGENCY  NAME  • AOORESV"  tUllatmtl  I 


i Control  I tn£  Oltlco) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  B »dSt.  yNlY  NUMBERS 


Project)  L 1^2622 AS  54fV 

'll,  AMtOKL  IUTS-. 

[\\  pp  II  J 
naalzsp.j 

IS.  SECURITY  CLASS,  (ot  thlfnport) 


UNCLASSIFIED 

la.  DECLASSIFICATION/ downgrading 
schedule 

NA 


I IS.  DISTRIBUTION  STATEMENT  (ol  Mia  Ktpori) 


{iv  siixi  • 


V M 4 


Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (at  tho  abstract  on  t trod  In  Bloc*  20.  It  dll  form  t I 


ft  tr  3 0 ' 

IS.  SUPPLEMENTARY  NOTES 


At 


1*.  KEY  WORDS  (Contlmn  on  nnni  It  n«c**M>7  and  Idnlliy  * y 6 lock  numbor) 

(U)  Gyroscopic  stability 
Dynamic  stability 
Equilibrium  roll  rate 
Resonance  roll  rate 

\Simple  particle  trajectory 


1 10.  A BA  TRAC 


mm  pmwmrm  1 1 1 


l HtmmUtr  bf  lf«dt  numtbmt) 


(U)  'The  report  describes  a computer  program  that  predicts  the  flight  trajectory,  spin 
decay,  and  gyroscopic  and  dynamic  stability  factors  for  spin*  or  fin-stabilized  cannon-launched 
projectiles  with  liquid  or  solid  payloads,  . 


DO  ij£Tn  1473  EDITION  OF  7 NOV  ••  IS  O DEPLETE 


y/y  j7<f 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  TNI*  PAGE 


smbNN4vj 

T 


PREFACE 


The  work  described  in  this  report  was  accomplished  as  support  effort  for  several 
projectile  research  and  development  programs  being  conducted  at  Chemical  Systems  Laboratory. 
The  work  was  started  in  September  1976  and  completed  in  September  1977.  The  computer 
program  is  general  in  nature  and  should  be  applicable  to  any  study  of  a projectile’s  trajectory. 

Work  was  authorized  under  Project  1L162622A554-1,  Lethal  Chemical  Weapons 

Technology. 

Reproduction  of  this  document  in  whole  or  in  part  is  prohibited  except  with 
permission  of  the  Director,  Chemical  Systems  Laboratory,  Attn:  DRDAR-CU-R,  Aberdeen 
Proving  Ground,  Maryland  21010;  however,  DDC  and  the  National  Technical  Information 
Service  are  authorized  to  reproduce  the  document  for  US  Government  purposes. 


SUMMARY 


This  report  describes  a computer  program  that  predicts  the  flight  trajectory,  spin 
decay,  and  gyroscopic  and  dynamic  stability  factors  for  spin-  or  fin-stabilized  cannon-launched 
projectiles  with  solid  or  liquid  payloads.  The  program  provides  a fast,  accurate,  and  efficient 
method  of  evaluating  projectile  trajectories.  A short  introduction  provides  the  mechanics  of 
trajectory  theory.  However,  no  effort  is  made  to  give  a complete  mathematical  analysis  of  the 
theory.  The  theory  may  be  found  in  the  literature  cited  (given  as  footnotes  in  the  report). 
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A COMPUTER  PROGRAM  FOR  ANALYZING  PROJECTILE  FLIGHT  TRAJECTORY 
I.  INTRODUCTION. 


The  purpose  of  this  report  is  to  present  a computer  program  that  assists  design 
personnel  in  analyzing  projectile  flight.  After  aerodynamic  coefficients  for  a proposed  munition 
shape  are  derived  from  wind  tunnel  tests,  a prediction  of  projectile  stability,  range,  and  accuracy 
must  be  made.  In  addition,  the  effect  on  projectile  flight  of  a liquid  payload  or  shift  in 
center-of-gravity  location  may  be  required.  This  program  can  be  used  as  a design  tool  in 
analyzing  projectile  flight. 


II.  THEORY. 

A.  Simple  Particle  Trajectory. 

As  in  any  trajectory,  simple  particle  trajectory  is  the  curve  traced  by  the  center  of 
gravity  of  a projectile  with  respect  to  a horizontal  surface  when  fired  at  a particular 
charge/elevation  combination.  However,  simple  particle  trajectory  uses  velocity  and  aerodynamic 
drag  to  calculate  the  horizontal  and  vertical  velocity  components. 

The  fundamental  equations  used  in  the  computer  program  are: 


< -drag 

- g sin  0 ^ At 

(1) 

v projectile  mass 

/ -g  cos  0 N 

*-(  v , 

)At 

(2) 

AX  - (V  cos  0)  At 

(3) 

AZ  = (V  sin  0)  At 

(4) 

where 

V = projectile  velocity 

0 = acute  angle  between  a horizontal  plane  and  the  tangent 
to  the  trajectory  at  the  projectile  center  of  mass 


g = gravitational  constant 
t = time 

X = horizontal  distance 
Z = height  above  horizontal  plane 
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V (a  FACTOR) 


When  a projectile  is  fired  from  a gun,  the  liquid,  initially,  is  not  rotating  with  the 
projectile.  During  the  down-range  flight,  the  liquid  acquires  angular  momentum  from  the 
projectile  and  the  liquid  spin  rate  approaches  that  of  the  projectile.  This  process  of  the  liquid 
acquiring  angular  momentum  is  referred  to  as  the  spinup  process.  This  process,  whereby  the 
liquid  acquires  angular  momentum  from  an  initial  condition  of  zero  rotation,  was  explained  by 
Wedemeyer.*  Wedemeyer  showed  that  for  the  case  of  a completely  filled  cylinder,  viscous  regions, 
known  as  Ekman  layers,  near  the  end  walls  act  as  centrifugal  fans  sucking  non-rotating  liquid 
from  the  in  viscid  core  and  throwing  the  liquid  outwards. 

To  incorporate  liquid  effects  (spinup  process)  on  rigid  body  spin  decay,  the 
difference  between  rigid  body  rotation  and  zero  liquid  rotation  at  the  muzzle  is  calculated.  The 
theoretical  spin  decay  assumes  a percentage  of  liquid  rotating  with  the  projectile  at  muzzle  exit. 
The  particular  decay  curve  (see  the  figure)  for  the  percentage  of  liquid  rotating  with  the 
projectile  at  muzzle  exit  is: 

y = a1  (5) 

where 

a = percentage  of  liquid  rotating  with  the  projectile  at  muzzle  exit, 
t = time  after  launch 
y = spin  decay  correction  factor. 

The  decay  curve  asymptotically  approaches  the  t-axis.  The  value  of  y (equation  5)  i:  multiplied 
by  the  difference  between  rigid  body  payload  rotation  and  zero  payload  rotation  at  muzzle  exit 
to  determine  projectile  spin  decay  due  to  liquid  effects. 


6 6 7 8 9.  10  11  12  13  14  16  16 

t (SEC) 

Figure.  Spin  Decay  Curves  for  Liquid-Filled 
Spin-Stabilized  Projectiles 

* Wedemeyer,  F.  H.  BRL  Report  1252.  The  Unsteady  Flow  Within  a Spinning  Cylinder.  Aberdeen  Proving 
Ground,  MD  21005.  1963.  10 
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C.  Roll  Rate  For  Fin-Stabilized  Munitions.* 


Unlike  spin-stabilized  projectiles  whose  spin  is  generated  by  the  gun  tube  rifling, 
fin-stabilized  munitions  generate  spin  by  the  air  flowing  over  the  canted  fin  blades.  The  spin 
equilibrium  roll  rate,  pe,  for  fin  munitions  is: 


I 


where 


Pe 


"Cls 

C 


VSfc 

lpd 


(6) 


Cls  = roll  moment  coefficient  due  to  fin  cant  (at  zero  spin) 

Cjp  = roll  damping-moment  coefficient 
Sfc  = fin  cant  angle,  radians 
V = velocity,  ft/sec 

d = maximum  projectile  body  diameter,  ft 

Theoretically,  resonance  instabilities  can  occur  in  spin-stabilized  projectiles  due  to 
the  coincidence  of  spin  and  yaw  frequencies,  however,  this  phenomenon  is  more  likely  to  occur 

with  fin-stabilized  projectiles.  The  spin  is  most  likely  to  coincide  with  nutational  frequency  0 j, 
given  by: 


where 


» - (1  + (h/sr)0-5 
2 B 


Sg  = gyroscopic  stability  factor 
pSd3  _ 

M = 21  CMa 

Z1y 

CMa  = static  moment  coefficient,  per  radian 
d = maximum  body  diameter,  ft 
p = air  density,  slug/ft^ 

S * projectile  cross-sectional  area,  ft^ 


(7) 


i 


■ 11  — — “ 

• Engineering  Design  Handbook  AMCP  706-242.  Design  for  Control  of  Projectile  Flight  Characteristics. 
Headquarters,  US  Army  Materiel  Command,  Washington,  D.C.  1966. 


Ix  . pd  i 
P = - and  <*>,  |=  “ 

ly  V j 

u = static  moment  factor  (per  radian) 
p = spin  in  radian  per  caliber  of  travel 
Ix  = axial  moment  of  inertia,  slugs  - ft^ 

Iy  = transverse  moment  of  inertia,  slugs  - ft^ 
This  reduces  to  a resonance  roll  rate,  p_,  represented  as: 


(rrj 


It  has  been  suggested  that  the  roll  rate  at  exit  for  a fin  munition  fired  from  a rifled 
tube  should  be  three  times  resonance  roll  rate.  A fin  munition  fired  from  a smooth-bore  tube 
exits  with  zero  spin.  Theory  states  that  munition  roll  rate  will  pass  through  resonance  roll  rate 
prior  to  equilibrium  roll  rate.  The  greater  the  equilibrium  roll  rate/resonance  roll  rate  ratio  the 
shorter  rime  the  projectile  is  in  the  vicinity  of  resonance  roll  rate,  thus  temporary  growth  in  yaw 
due  to  resonance  will  be  insignificant. 

III.  THE  PROGRAM. 

The  program  given  in  Fortran  symbols  on  pages  13  through  20  calculates  the  flight 
trajectories  of  spin-  and  fin-stabilized  projectiles  with  liquid  or  solid  payloads.  The  input 
parameters  (i.e.,  launch  conditions,  aerodynamic  coefficients,  projectile  physical  characteristics) 
required  are  listed  at  the  beginning  of  the  program.  The  program  consists  of  three  parts  (l)the 
main  program  calculates  trajectory  and  stability,  (2)  a subroutine  calculates  the  liquid  effects  on 
projectile  spin,  and  (3)  a subroutine  calculates  equilibrium  and  resonance  roll  rates  for  a 
fin-stabilized  munition.  Sample  outputs  for  solid  and  liquid  spin-stabilized  projectiles  are  shown 
in  tables  1 and  2.  Table  3 shows  an  output  for  fin-stabilized  projectiles. 
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201 
2 07 
2 P3 
205 
7 0S 


inrjrjn 

50  0 IFIJ.LT.KI  60  TO  69 
VRITEI6.5SI 

55  TORN  AT  (IX*  20H  ERROR  STATEMENT  5001 
STOP 

11  1 CO=CDT  f JI 
CMrCMAI JI 
SA=5A(5(  J1 
oa=damiji 
CSL=CL5IJI 
aop;cLPoi  ji 
a - CLA (JI 

IF  I OPT. GT. 01  SO  TO  5 0 
GO  TO  81 

11  2 C0=C0T(J-1  I♦(C0T(J1-C0T<  J-l  II  / 1 VTI  J l-V  Tl  J-|  » »•  IV -V  T(  J-l  I I 
CHtCMA I J-l I +<  CMAI J) -CMAt J-I II /< VT( JI-VTI  J-1 I I* (V-VTI J-l I 1 
SA=9AG( J-l ! ♦( SAGI JI -SAG! J-l 1 1 /( VTI Jl-VTI J-l I I* IV -VTI  J-l I > 

OA  :0  AM  I J-l I ♦( 0 AMI JI -QAMt J- 1 1 1 / 1 V Tl J I -VT I J- 1 I 1 • I V -V T I J -1  I I 
CSL  = CLS(  J-l )♦ I CL5I JI-CL5I J-ll 1/ I VT I Jl-VTI J-l ||  •(  V-VTI  J-ll  l 
CL  0°:  CLP  01  J-l  !♦  (CL»D(  Jl-CLPOI  J-  1 11/  I V’  I JI -VT  I J-l  1 1 * I V- VT  I j-  1 > I 
CL:riA I J-l I ♦ICLAI JI -CL  A I J-l I 1/1 VTI JI-VTI J-l t I* (V-VTI J-l I I 
IF  (OPT.GT.OI  GO  TO  90 

jl  NCNT  1 :0 
NCNT2=0 

C VELOCITY  '•OPPONENTS 
VX:V*CC5ITHETA) 

V2  = V*S INITHET  Al 
FACTOR  : .S*CiLN*cn*S 
GO  TP  163 

50  FACTOR  = . 5 *VOE  N*  C 0*5 
C DRAG  C OM°Q  N£ NT  5 

159  DR  AG  = FACTCR*V *V 

OX:ORAG*COS  tTHfTAI 
GZ:0RAG*S1M  THETAI 
C ACCELERATION 

AXt-DRAG/WT  • CO S I THETAI 
AZ=-OR AG/ «T  *5  IN  I THETA  1-32. 17 
THE 4CC  = 4 TAN2I  A7 . AX1 
A:909T I AX*AX* A2 • A 2 1 
IF  I 2 t 50.50*45 

55  THTACCrTHE  ACC*1  80.2  3.  15  153 
THrT  A:  THET  A*1  «0./  $.14  15  3 
C PMNT  OUT  INTERVAL 

IF  (NCNT  l.NE.  NCNT2I  GO  TO  ?1 
IFIJD.EO.il  GO  TO  360 

735  FORMAT (4F1 5.4  1 

SPR=VO/ (01 A*TVI5T0> 

5P  INrSPR 

IFIJA.EG.OI  GO  TC  360 
SPINL:V0/(0IA*TVI5TLI 
SPTNZSPINL 

360  IFIJC.GT.OI  GO  TO  300 

VHITEI6.2I  T*THETA.Xl.2l.V.A.THtACC.09AG.C0«Vr«r»'. 

1SPIN.SPRC1 

2 FORMAT  II  M,  2X.5F  3.2*  2F  3.5*  IF  j.  2*  2*7.  31 
GO  t0  J10 

JO  0 mfit  I TE  I 6 . 4 t T*THETA.X1*71.V,A  , THE  AC  C » OR  AG  . C 0.  VO  r V . S <>I  N 
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fv  rj  'g  (g  <\.  fj  fg  w fw  \j  i\j  l j ‘\j  / v >\*  f\J  f\/  r~j  t\J  /V  ’ J > « r*»  r\j  r>j  i\>  ;\>  i j i * K/  K/  fw 


TO^ 
2 07' 
208 
2 09 
10 
11 
12 
1 3 
14 
19 
16 
17 
13 
19 
iO 
21 
22 
23 
^4 

25 

26 
27 
23 
29 
70 

31 

32 

33 

34 
.>5 

36 

37 
39 

39 

40 

41 
*•2 

43 

44 

45 

46 

47 
2 48 
2 49 
2 50 
251 
25’ 

253 

254 
2 55 
2 56 
2 57 
2 58 

259 

260 
261 
26? 


4 ?6PM*r  rfTT8F9.2  »2F9.5»  1F9.2I 
310  NCNT1  = NCNT1  *INC 
2 1 NCNT 2=NCNT  2*1 
VXSAVE=VX 

ves«vc=v2 

VSAVE=SQ  RT  (VXSAVE*VXSAVE*V?SAVE*VZSAVE) 
THETAS=ATAN2fV2 SAVE. VI  SAVE! 

TSAVF=T 

XSAVFrX 

2SAVF=2 

OXSAVE=OX 

N5AVE=DZ 

AXS A VE- A X 

AiSAVEtAZ 

c increhenting  trajectory 

DELVX=AX*OELT 
DEL  V2  = A?*0ELT 
DELX:DELT*tVX*.  5«DELVX» 

DFL?:QELT*t  VZ».  5«0ELV2> 

VX: VX+DELVX 
V2=VZ ♦DELV? 

T=T ♦OELT 

THE T*:AT AN21VZ.VX 1 

V=6QRTtVX*VX»VZ*V2l 

X=X*OELX 

X0I5T  rX-XS  AYE 

Xl=X/3.281 

2 = 7 ♦ CELZ 

ZQTST :?-ZS  AVF 

YOFN  z »LXPH-2HsL-21 /?2not.U»l  *DEN3L 
IF  (VOFN.LE  .DEN)  SO  TO  <*3 
VOr N : OEN 
4 7 21-2/3.281 
C SPIN  CALCULATION 

FAC =VOEN*. 735*01  A *01  A/I /.0«VT» 

16  ( JO  . EO  . 0 1 60  TO  214 
CALL  FTN 

21 4 OISTrSORTt XQIST»XOIST»?DIST*ZOIST> 

RG  A =V  8*  D IA*  *2  ./  AXN  0 
RGT  rw  a*  D I A*  *2  ./ T9  A N 
ARATlo:COS«THET  A)  / ( C051  THETAS  »» 
a°=-CD/»3.0*RGA» 

C06=RGA*CLP»C0 

SRPC  = SRPCO» ARATIO-I EXP»FAC*OIST*COE  » 1 
TWI6T0=I6.28318*V0J/IS9PC*VJ 
SRpC0=$9  pC 

IFtJA.E9.0l  GO  TO  370 

IFITI.Gf.TMAX  1 GO  TO  37  0 

0SRPC=S9PC* f AXM/AXM01 

SP  1=SRPC-0SRPC 

CALL  DECAY 

SP2=W*Spl 

SPL=SP?»DSRPC 

TW ISTL = <6.283 18 «VOt /IS’L»V» 

IF  t TI.PC.THAX » GO  TO  37 1 
GO  TO  37U 

3’  1 IF  tFCPT.EG.  2*1  GO  TO  370 


J 
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f 

■j 


263  • 

WiSTd=TNlSTL 

264  • 

SR*CO  = SPL 

2 65  * 

JA=0 

266* 

K0PTi2 

267* 

C 

YAM 

CALCULAT  ION 

268* 

370 

AOEN=VOEN/32. 17 

269* 

SPRC=SRPC*(  U/DIAl *V 

27(3* 

YAWl=AXMS*SPRC*  32. 17*C0SI  THETA! 

271  * 

YAW2=.5*ADEN*S*DIA*CM*V**3. 

272* 

YAW3=t VAW1/VAM2 1*180. /3. 14159 

2 7 J * 

YAW4:(SA*SPRC*0IA)/(CM*V) 

274  • 

VAW5=(YAU3*YAW4 1*1 80. /3. 14159 

275* 

SP  R C 1 = CM  PL  X ( Y AW  3 • Y A MS  1 

2 76  * 

IF  1 OPT. G T.O  1 GO  TO  35 

277  * 

C 

STATIC  MOMENT  FACTOR.  L«-FT/RAOIAN 

278  « 

U=  .3926 9* ADEN *01 A** 3. *V**2.  *CM 

279* 

IF  < JB.GT.Ol  GO  TO  3 5 

?80  * 

IF(JA.GT.O)  GO  TO  36 

281  * 

C 

SO  LI  0 FTLLEO  PROJECTILE 

282  • 

C 

GYROSCOPIC  STABILITY  FACTOR 

2 8 3 * 

SG-(AXM$**2.*SPRC** 2. ) / (4 . *TRAN  *U) 

284  * 

C 

DYNAMIC  stability  estimate 

2 e 5 • 

S01=2.*<CL*RGA*SA> 

286  * 

S02=0L-CD+1 -rgt*da> 

287  • 

GO  TO  37 

?86  • 

C 

LI  QU  10  FILLED  PROJECTILE 

28  9 • 

36 

SG1=(AXMS**2. *S°R  C*  *2 .) / C 4.*TRAN*Ul 

290* 

SG2=AXM9/ AXMS 

291  • 

SG=SG1*SC2**2 . 

292* 

»GAP=(WP*DI»**2 .l/AXMR 

2 9 i • 

RGTR=«WW*OIA*  *2  .»  / TRAN 

2 94  * 

S01=2.*«CL*RGAR*SA) 

2 95* 

S02:CL-C0-(RGTR*CA» 

2 96  * 

3 7 

SOOrSO 1/S02 

297  • 

SG1=1./SG 

298  « 

SD=SOO*<  2.-S00) 

299* 

35 

0PT  = 1 

300* 

GO  TO  16 

301  * 

C 

OETFRMINE  POINT  OF  IMPACT 

3 02  • 

50 

«ATIO=-ZSAVE/OFL? 

303  * 

XF  INAL  =R  AT  1 0*  OE LX*  X S AVF 

304  * 

XMFINL=XFINAL/3  .2B1 

305* 

ZFINAL=P  ATIO*QELZ*ZSAVF 

3 06  • 

2*FINL  = 7FINAL/3  .281 

307* 

TFTN«L=3ATI0*  « T-TS A V E >*T SA Vc 

''08* 

VXFNL=RATIO*DELVX*VXSAVE 

3 09  * 

VZFNL=RATIO*OELVZ*VZSAVE 

310* 

VFINAL=S«T»VXFNL*¥XFNL*V2  FNL*V7FNL» 

3 11* 

THEFNL=ATAN21  VZFNL.VXFNL) 

312* 

AXFNL=RATIO*( AX -AX  SAVE) *AX SAVE 

313* 

AZFNL: RATI  0*1 AZ -AZSAVEl *AZSAVE 

3 14* 

AFINAL=SORT  C AYFNl *AXFNL*AZFNL*AZFNL  » 

315* 

THACFL=A  TAN2I AZFNL. AXFNL1 

316* 

DXFNL=RATIO*l 0X-0XSAVE1 *OXSAVE 

317* 

KFNL  = RATIO*«OZ-DZSAVE1*OZSAVE 

3 18* 

OF  INAL  =SO«T(OXFNL*OXFNL*OZFNL*OZFNL» 

319* 

THACFL=THACFL*1  8U./3.  14159 

im 
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320  • 

THEFNI.rTHeFNL.HO. / J. IMS's 

121  * 

VRITEI6*  31  TFINAL.THEFNL.XMFINL.2HFINL.V  FINAL* AF INAL. 

12?* 

ithacfl.ofinal 

3 21* 

1 

F0RNAT11X.15H POINT  OF  IMPA CT . /.  IX .8 F9 .2 » 

5 2*4  * 

C 

CDMPARF  SG  TO  SOO 

32S. 

If! JO. E8.01G0  TO  J79 

3 * 

VRTTFI6.381  IPPATIO 

2 2 7 * 

1o  1 

FORMAT!/. IX. 42H0ATI0  OF  E<JUI  LIBRIUM  RCLL/RESONANCE  ROLL  = .F 

5 2 f • 

179 

IF  1 HOP  T. NE . 2)  GO  TO  471 

M )• 

JA  r 1 

* ’0* 

47  1 

IF  1 JS . GT. 01  GO  TO  49S 

3 31  * 

IFIJA.ra.l 1 GO  TO  10G 

5 32  * 

HRITEI6.  1?1  1 

’33* 

12  1 

F0PMATf/.lX.2  3HS0LID  FILLED  PROJECTILE! 

3 ju  • 

GO  TO  122 

3 33  * 

12  6 

NRITEI6.11U 

3 3b* 

11  3 

FORMAT  ( / . 1 X .2  4HLI3UID  F I LL  rP  PROJECT  TLFI 

3 37  * 

12  2 

Nfl  I T£  1 6 » 101  1 SG 

11*  * 

12  1 

FO  PM  AT  1/. IX. 3 OH  GYROSCOPIC  LABILITY  FACTOR  = .nn.*l 

l 2 * * 

JR ITE ( 6 « 102  1 SOO 

mn  • 

10  2 

FORMAT!/. IX. 27HQVNAMIC  STABILITY  FACTOR  r .F13.4J 

3 4 1 * 

IF  I SGI . GE .SOI  GO  TO  2 10 

3 4 ? » 

WR  T TE ! 6 » 20  5 ) SGI. SO 

143* 

.0  4 

FORMAT  1 /. IX . 4 9H PROJECT  1LE  IS  STABLE  SINCE  S l./'5  < SUCI2.0 

144* 

1F10.4.3H  < .F7.41 

344* 

GO  TO  4«5 

3 46  * 

21  0 

•RITEI6. 220 1 SGI. SC 

347  * 

22  0 

FORMAT!/. IX. SlHPRO JECT ILt  IS  UNSTABLE  SINCE  5 1. /So  > SD0!2 

34c » 

1 .F in. 4. 3 H > .F7.4 I 

3 4 3 * 

4,6 

CONTINUE 

3a1* 

4 

53  a 

END 

1 • 

c 

THIS 

SUBROUTINE  CALCULATES  THE  PROJECT  ILF  S»IN 

? • 

c 

OE  CA  Y ASSUMING  VARIOUS  LIQ  SPINUP  IN  THE  TUBE*9* 

3 * 

c 

4 IS 

THE  PERCENT  OECAY  AT  TIKE.TI.  FOR  A GIVEN  R 

4 • 

SUBROUTINE  OECAY 

s • 

COMMON/ 8LK 4 /B .N.TT.TMAX 

6 * 

IF!8.GT.0.0.AN0.3.LT.l.n>  GO  TO  10 

7 • 

NRITEI6.5) 

8 * 

5 

FORMA T 1 1 X . 2 1H  ERROR  PERCENT  SPINUPI 

9 * 

STOP 

in* 

1 0 

NrB**TI 

n * 

TlrTI*.02 

12* 

IF1TI.LT.TMAX 1 GO  TO  20 

15* 

TMAXrTI 

14. 

GO  TO  70 

IS* 

2 0 

IF  IN. GT. 0.001  1 GO  TO  70 

16* 

TMAXrTI 

17* 

70 

RETURN 

19* 

ENO 

19 


1 * 

C 

subroutine  calculates  spin  pop  a fin  munition 

2 * 

c 

PE  IS  THE  E0UILI8RIUH  ROLL  9 AT f . 

3* 

c 

PR  I 

S THE  RESONANCE  ROLL  RATE 

4 • 

SUBROUTINE  FIN 

5* 

COMNON/BLKl /C SL .SF C . CL  0° . 0 PT . P. T» 

a • 

COMMON/ RLM  2/DEN*  V*S «DIA*D£LT 

7 * 

COMMON/ 8 LK 3/CM. SPIN* PR  AT  10 . AXMS * TRAN 

o * 

REAL  MU 

9 * 

OOUBLC  PRECISION  SPI.E 

n* 

E=2.7182G**T 

n * 

Oi=-CSL/CLDP 

12  • 

D2  = <V*SFCI/QIA 

13  * 

SPIN=01*D2 

14  * 

IFfT.EQ.S.)  GO  TO  50 

1 5 • 

IFIP.EO.SPIN)  GO  TO  SO 

t n * 

IF  1 *»LT.  SPIN)  GO  TO  40 

17  • 

SPI=SPIN*(P-SPINI*{ l./E» 

1*  * 

GO  TO  45 

n * 

40 

SP I = »♦(!.-! 1. /E 1»  *(SPIN- PI 

2C* 

45 

SPIN:SPI 

cl  • 

50 

MU=.5»DEN*V**2. *S*DIA*CM 

22  • 

IF  fOPT.GT.O  . 1 GO  TO  30 

2 3 » 

IF  ( CM.LE.Q.  0)  GO  TO  10 

i 4 » 

UR 1TE (6*  20 1 

25  • 

20 

FORMAT! /.  IX  .45H  SUBROUTINE  7 f OUIRE  S NEGATIVE  STATIC  MOM  COEF» 

?'  • 

STOP 

27* 

1 0 

PR  -1  SORT (-MU/ (TRAN- AX MSI  1 ) / 1 2 .*  3.  1 4 15  ) 

23  * 

PRATIO=SPIN/PP 

29* 

30 

RETURN 

37  * 
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Table  2.  Output  for  Spin-Stabilized  Liquid-Filled  Projectile 
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IV.  CONCLUSIONS. 

The  program  analyzes  projectile  flight  as  shown  by  the  sample  outputs.  These 
predictions  for  spin-  and  fin-stabilized  projectiles  will  assist  design  personnel  in  evaluating  a 
projectile  flight.  However,  it  should  be  noted  that  since  the  program  does  not  take  into  account 
meteorological  conditions,  other  than  air  density,  and  that  it  adapts  simple  particle  trajectory 
theory,  the  output  provides  a good  trajectory  estimate  for  a projectile  fired  under  normal 
conditions. 
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